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^ ■ Abstract 

a^ , 

Q^ , It has been proposed by Cook {Phys. Scr. T 21, 49 (1988)) to use a short probe laser 

pulse for state measurements of two-level systems. In previous work we have investigated 
to what extent this proposal fulfills the projection postulate if ideal photon detectors are 
considered. For detectors with overall efficiency less than 1 complications arise for single 
systems, and for this case we present a simple criterion for a laser pulse to act as a state 
measurement and to cause an almost complete state reduction. 
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O ■ 1 Introduction 

• The outcome of measurements is often described in an idealized way by the projection postulate 

0> . with its state vector reduction. The current form of the projection postulate is due to Liiders 

r~| [ |||] and for degenerate eigenvalues it differs from that of von Neumann . For non-degenerate 

Qh' eigenvalues it is implicitly contained in Dirac's book For the simple case of a two-level 

system in a state and for a measurement of the energy the projection postulate states that, 
Q i with probability KljV')!^, one finds Ei and the state is changed to correspondingly for £^2- 

In this simple case a measurement of the energy can also be regarded as a measurement of the 
eigenstates |1) and |2). 

Cook Q proposed to perform a measurement of the states |1) and |2) of such a two-level 
^ I system by means of a short laser pulse which drives the transition between level 1 and an 

' additional, rapidly decaying, level 3 (see Fig. 1). For a single system, initially in a superposition 

of |1) and |2), the short laser pulse may, or may not, lead to the emission of fluorescence photons. 
In the former case the system was taken to be in |1) after the end of the pulse (and after a short 
transient decay time for level 3), since right after each photon emission it is in the ground state. 
In the latter case the system was taken to be in |2). This was used in an experiment on the 
quantum Zeno effect 

In two previous papers |^] we have investigated in how far this proposal realizes a 
measurement to which the projection postulate can be applied. Related work appeared in 
Refs. ^, |ll], 12, 13]. It was shown in Ref. |^] that, under the condition of no photon emis- 
sion during a pulse, the reduced (or conditional) time development |]lj, F7, 15] carries the 
state of the system ever closer to ]2). The quality of this reduction to ]2) can be expressed by a 
single combination of the pulse parameters, namely Ne{Tp; ]1)), the average number of photons 
emitted during a pulse of length Tp when one starts from )1) as initial state. can be calculated 
from the fluorescence rate of an ensemble of such systems. It was shown in Ref. that for A'^e 
as low as 5 or 10 one has an excellent reduction to 12). 
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In quantum mechanics it is in principle not legitimate to speak about emission of a photon 
unless it is observed (detected). However, it was pointed out in Ref. [14| and Ref. Ill] that in 



questions of photon statistics it makes no difference theoretically whether or not the photons 
are monitored over all space at rapidly repeated time intervals. Emitted photons are then those 
which would be detected by ideal detectors in the whole space. 

The discussion in Refs. ^ is based on such ideal detectors. With non-ideal detectors 
problems arise in the reduction of the state of a single system when no fluorescence photons are 
detected during a pulse, because this null detection may be due either to no emission at all or 
to a failure of the detector to record emitted photons. In the former case the system is in |2) at 
the end of the pulse, while in the latter case it is in Hence the system has to be described by 
a density matrix with a |2)(2| and a |1)(1| contribution. For a good reduction to |2) the latter, 
the "non-reduced part" should be small, and the aim of this paper is to give a simple criterion 
to ensure this. 

A real detector may miss a photon because of low quantum efficiency or because it operates 
only in a smaller solid angleQ. We define the overall efficiency r/' as the ratio of detected vs. 
emitted photons. 

For a system initially in the ground state, we denote by A'^d(rp;|l)) = 77'A'c(rp; [1)) the 
average number of fluorescence photons detected during a pulse. Experimentally this is easy 
to determine from the observed fluorescence of an ensemble ("gas") of systems, and one can 
adjust the pulse parameters to realize any value of for a given detector. In a simplifled 
derivation we will show in Section 3 of this paper that for initial state IV") the non-reduced part 
is approximately 

1(11^)1' e-^d n) 



for small overall efficiency r]'. However, this simple derivation gives no insight on how good the 
approximation is and how small rj' has to be. Therefore we show in the precise treatment of 
Section 4 that the non-reduced part is smaller than 

Kl|^)|2(l+ ^ V)e-^<^ (2) 

for r]' < 1/3 and > 1, provided is sufficiently large, e.g. 5, to ensure reduction in the case 
of ideal detectors. This is a strict bound from above; a strict bound from below is given at the 
end of Section 4 in Eq. ([57|). If the initial state is a density matrix p instead of a pure state I'ip), 
then Kll'i/;)!^ is replaced by pn. We derive Eq. (§) by Laplace transforms. Similar techniques 
have been used in Ref. l2^ 



2 Basics and previous results 

We consider a V system as in Fig. 1, with level 2 (meta-)stable and level 3 rapidly decaying 
(Einstein coefficient A). The laser pulse has Rabi frequency Q and is tuned to the 1-3 tran- 
sition. In the quantum jump approach (or Monte-Carlo wave function approach, or quantum 
trajectories) []l4| p^ , p^ , 17, 18, |2^] it is shown that under the condition of no photon emission 
the time development of a system is given by a reduced (or conditional) Hamiltonian. In our 
case it is given in an interaction picture by Q 

iJrcd = -i?iM = I {1^(|1)(3| + |3)(1|) - iA |3)(3|} . (3) 



^ Solid angles smaller than Atv have been considered in Ref. [ po| , where in one solid angle photons were 
counted, while in the other solid angle the light was spectrally analyzed. 
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The eigenvalue of M are A2 = (with eigenvector |2)) and 
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Ai,3= 4 (At - 41^2 j . (4) 
For initial state 

IV') =ai|l) + a2|2) , (5) 
the conditional time development is then 

e-'^-^/^|V')=a2|2) + ai j^^^^^^^^^f 11) + ^ ^ _ |3)| • (6) 

The probability for no photon emission until time r, denoted by -Po(t; is given by 

Po^(r;|V')) = ||e-'^'--/^|V')ir (7) 

From Eq. (^) one obtains 

Po^(r;jV')) = 1(21^^)12 + [(l|^)|2pe(r;|l)) . (8) 
The probability density for the first photon emission is 

wl{r-M)) = -^Po'(r;|^)) 

= \m\'wl{T-\l)) (9) 

since the 1(21-0)1^ term in Eq. (|8[) drops out. In particular, for an atom with initial state |1) one 
has 

and 

If Tp, the duration of the pulse, satisfies 

Tp » max{l/A ; ^/Sl^} (12) 

then for r = Tp the exponentials in Eq. (^) can be neglected and the state is proportional to |2), 
with probability 

p^,{T^-M)) = \m? , (13) 

just as with the projection postulate. In this way one can understand why a system without 
emissions during a pulse ends up in the state |2). 

With the complementary probability 1— 1(21-0) p = |(l|'0)p an atom will emit photons during 
a pulse and will be in a superposition of |1) and |3) at its end, at time Tp. The |3) component 
will then either lead to the emission of an additional photon or the reduced time development 
(now with = 0) leads to its vanishing (without emission). In both cases the state will be in 

after a transient time rtr whose duration is several multiples of ^4"^. 

Analogously, for a system with no emission until Tp + Ttr the tiny |3) component, still present 
at time Tp, has vanished by then and the remaining tiny |1) component of exp{— ii/j-ed'^p/^llV') 
in Eq. (^) is the part which is not reduced to |2). This can be estimated more precisely in terms 
of A'c(Tp; |1)). The latter is obtained by integrating the photon rate Ap33{t) where pss^r) is 
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given by the Bloch equations for the driven 1-3 two-level system with pii(O) = 1 [^]. This 
gives 

'v.(-; ID) = [Ar - + ■ j (14) 

where 

^1/2= ^ (SA T _ i6j^2) . (15) 

For A'^c > 1 the contribution of the exponential terms is negligible. (For weak driving the time 
for the first emission is long which makes exp(— /ijr) small; for strong driving the factor in front 
of the last term is small.) We can and will therefore use, for A'^c > 1, 

q2 / 3^2 \ 

With this one can express Tp by Ne{Tp; |1)) (and the remaining parameters) and can insert it for 
r on the r.h.s. of Eq. (^). A graphical evaluation then yields as upper bound for the non-reduced 
part 

|(lIexp{-ii7,edTpA}lV)P < imi' e-^^(-p;|i)) . (17) 

Thus in the case of no photon emission during a pulse the reduction of {ip) to |2) is already 
excellent if the parameters of the pulse are such that iVe(rp; |1)) > 5, i.e. if the average number 
of photons emitted during the pulse by an atom initially in the ground state is 5 or more. The 
estimate in Eq. (p^) can be improved substantially for strong driving 



3 Simplified derivation 

We now incorporate a non-ideal detector of overall efficiency ij' < 1. With initial state \^) as 
in Eq. (^) we denote by Pq{t; lip)) the probability of no detection of fluorescence photons until 
time r for a single system. With probability P(j{Tp; \Tp)) no photon is emitted during a pulse, 
and we assume in the following that for ideal detectors the non-reduced part is negligible, i.e. 
the corresponding state is |2). With probability Pq — Pq photons are emitted, but missed by 
the non-ideal detector, and in this case the state is |1), after a transient decay time Ttr for the 
3 component. Hence a system, for which no fluorescence photons are detected, is described by 
the density matrix 

/,0(rp;r?') = Po'(rp;|V)) \2){2\ + [pi{Tp;m - P^irp;m) |1)(1| (18) 

where the normalization is such that the trace is Pq, i.e. the probability of no detection. We 
call the |1)(1| component the non-reduced part, and the smaller this is the better a reduction 
as in the projection postulate is realized. If one makes no photon measurements, i.e. r]' = 
and Pd = 1, then p^{Tp; rj') is just the diagonal density matrix of the Bloch equations which was 
discussed in Refs. p^] . 

We will show further below (Section 4.1) that for Pq an equation analogous to Eq. (^) holds, 

pi{T-m = \m? + \m?p^{rm . m 

The density matrix in ( [T^ ) then becomes, with Eq. (^) for Pq, 

pO(Tp;^') = K2!V')P |2)(2| + |(l|^)ppd(^p.|i)) |i)(i| . (20) 
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To determine the non-reduced part one therefore only has to calculate Pq{tp; |1)), the prob- 
ability to detect no photon for initial state |1). This will be done in detail in Section 4. Here 
we use a simplifying assumption, namely that for initial state |1) the fluorescence photons are 
emitted with a constant probability rate / during the pulse, 

/ = iVe(Tp;|l))/Tp . (21) 

For the detector the recorded rate is then 

=r?7 = iVd(rp;[l))/rp . (22) 

Similar to radioactive decay Eq. ( p2| ) now implies 

Po"(rp;|l)) = e-V-p 

= e-NArp;\i)) . (23) 

Instead of a constant rate one can also assume the emission process to be Markovian. This 
assumption leads to the same result (cf., e.g., Ref. [^, p. 231). 

We expect this to be a good approximation for small rj' due to the following. For small rj' most 
emitted photons go undetected, and to have a given average number N^^, = 5 say, one needs 
a large A'e- The detected number of photons in this time interval has then an approximately 
Poissonian distribution, in accordance with Eq. (p^). However, this argument gives no insight 
how small r]' has to be or how large the error is. To answer this we need the more detailed 
analysis of Section 4. 

From Eqs. (|23| ) and (^) one obtains Eq. (||). Thus, with the above simplifying assumption 
of constant detection rate the non-reduced part decreases exponentially with N^^Tp] |1)). For a 
given detector one can adjust the laser parameters accordingly to come as close to the projection 
postulate as desired. 



4 General derivation 

4.1 Solution by Laplace transforms 

With the results of Sections 2 and 3 the incorporation of detectors with overall efficiency rj' < 1 
is a purely statistical problem. After each emission the system (atom) is in its ground state, with 
all memory lost, and therefore the sequence of photon emissions is a renewal process [25|, as is 



the sequence of photon detections. This renewal process (although derived from the Markovian 
quantum trajectories of state vectors) is not Markovian since it would then be a Poisson process 



1 25, Hence the simplified treatment of Section 3 is only an approximation. 

In view of Eq. ( pO|) we have to determine Pq{tp; |1)) in order to estimate the quality of the 
reduction. For a system initially in {ip) the probability to have exactly N emissions occurring 
in (ri, Ti -|- dri), (tat, ttv + dr^) in the interval (0, r) is 

wlin; mdnwHT2 - ri; |l))dT2 • • • <(rjv - rjv-i; |l))drjvPo'(^ " ^n; |1)) (24) 

and hence the probability to have exactly N emissions in (0,r), denoted by P^{t; \tl')), is for 
> 1 obtained by integrating over the domain < ti < r2 < ... < TN <T. With Eq. (I) for 
Wi{ti; IV')) this gives 

PNir;m = r dTNP^{T-TN;\l)) r dTN-lwUTN-TN-i;]!)) 

Jo Jo 
'''dri^(r2-ri;|l))u;^(ri;|l)) . (25) 
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For N emitted photons the probabihty for all photons to go undetected is (1 — r/')^, and therefore 
the probability of detecting no photon at all in (0, r) for initial state is 

CO 

Po'(r; IV^)) = PSir; |^)) + E " v'fPl^ir; m • (26) 

N=l 

Eq. (19) immediately follows from this when using Eqs. ( p5| ) and (^). 

We now use Laplace transform, denoted by " . The convolution theorem for Laplace trans- 
forms together with Eqs. (|2^) and (p^) for |^) = |1) gives 



N=0 



(27) 



l-(l-r?Ou'f(p;|l)) ■ 

pQ and wf are easily obtained from Eqs. (|lO| ) and (lTl| ) and yield as a quotient of the form 

P^felD)^^ (28) 



where 



/.'(P) = Iap^+ ^{A^ + 2n^)p+ ^r^'An^ (29) 

5(P) = /+ 1(^2 ^21^2)= /o(p)/p. (30) 



The same expression was derived by somewhat different methods in Ref. |22| and the zeros of 
the denominator were found for the special case Q = A/ 2. 

Let pi, p2, P3 be the zeros of fr^'{p). Either all of them are real and negative or one is 
negative and the other two are complex conjugates of each other, with negative real parts. We 
choose pi to be the real zero closest to zero and | Rep2| < | Rep3|. By partial fractions one then 
obtains, for distinct zeros, 

^'o'(t;|i)) = E ^^^'^ 



1 {pi -Pi+i){Pi -Pi+2) 

= ciePi^ + C2eP2^ + cge^'^^ . (31) 

where we use Pi+3 = Pi- If two zeros coincide one has to take limits. An alternative form is 
obtained by noting that 

fr,' {p)/ {P - Pi) = {p- Pi+i){p - Pi+2) ■ (32) 
Letting p pi, the l.h.s. becomes f^i{pi) = foiPi), and thus 

Poir;\l)) = El^e^-. (33) 

i^l J0\Pi) 



Furthermore, since 



and f'q'ipi) = 0, one has from Eq. 



fo{p) = fn'ip)-lv'An^ (34) 



giPi) = -\rlAV?lpi . (35) 
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Again using f.q'{pi) = to recast Pif^iPi) one than finds 



(36) 

We will next derive an approximate solution to Eq. (|33| ) and will then discuss its range of 
validity. 

4.2 Approximate form of 

With /xi^2 as in Eq. (|l5|) the zeros of frj'=o{p) are 

{Pl,P2,P3) = (0, -fJ'i, -/i2) • (37) 

For small rj' , pi can be found perturbatively by linearization (Newton's method). In particular 
one finds 

This can now be inserted into Eq. ( |3^ ) for ci, where in the denominator we neglect r/'^ when 
compared to 1 and rj' . Eq. ( p!^ ) for Nc{Tp; |1)) is used to recast e^^"^?. With N^^ = rj'Ne we then 
obtain 

Expanding in r]' , the linear terms cancel and only ■q''^ terms with coefficients less than 1 remain; 
these are neglected when compared to 1. 

Prom Eqs. ( p!5|) and ( |37| ) one sees that, for small t]', e^^'^p and e^^'^pare much smaller in 
absolute value than e^^'^p. Therefore the second and third summands in Eq. (|3l| ) are neglected. 
(Some care is needed for p2 close to ps; this is discussed in detail in the next subsection.) 

Thus we arrive at the approximate expression for small r/', 

Po^(rp;|l))=e-^^(-P'|i)) . (40) 

In view of the derivation ij' has to be small enough for Newton's method to give a reliable pi. 
In the next subsection the range of validity of the above expression will be discussed in detail. 



The reader not interested in all the details can proceed directly to the strict bound in Eq. (55). 



4.3 Validity range and precise estimates 

We now present upper and lower bounds on Pq, and thus on the non-reduced part. These will 
justify the preceding approximate solution in more detail. 

For this we will need more information on the zeros of the curve y = fr^'{p). This curve is 
obtained by moving the curve y = fo{p) upward by the amount ^ 7]'AQ'^ or, alternatively, by 
moving the abscissa downward (see Fig. 2). We note further that, for rj' = 1, the zeros are, in 
the given order, 

/ (-2Ai, -i A, -2A3) for Ai real 
I (-i A, -2Ai, -2A3) otherwise 

where Aj is given by Eq. (^). The point of inflexion of the curve is at pinfl = — ^A, which is also 
a zero of fi{p), and this fixes the abscissa for r/' = 1. Thus, from Fig. 2, pi lies between p-^^a and 
0. On the other hand, pi lies to the left of the tangent at p = (see Fig. 2), and so the r.h.s. of 
Eq. (|38|), i.e. Newton's method, provides an upper bound for pi. One also verifies that 



(Pl,P2,P3) 
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by checking that (p) changes sign between these values. 
As zeros of the polynomial the pj's satisfy 

3 

P1+P2+P3 = -^A (42) 

PiP2+P2P3+Pm = ^{A^ + 2n^). (43) 

If P2 and p3 are real they must lie to the left of pinfi, i.e. p2,3 < — ^ ^- If Pi is the only real zero 
then p3 = P2- Writing 

P2/3 = -a^ib , a,b>0 (44) 
one obtains from Eqs. ( |4^ ) and (^) 

3.1 

a = 4 ^ + 2 ^1 

= f p?+ ^Api. (46) 

lb 4 4 



Now we discuss cie^^'^p, with ci > given by Eq. ( p6|) for i = 1. We replace in the 
exponent and in the denominator of ci by the r.h.s. of the inequality ( |4l| ) and by the l.h.s. 
of (PI). Then one obtains instead of Eq. (^), the inequality 

for ?/ < ^. Evaluating this graphically one finds, for arbitrary A and Q and for r]' < 1/3 

cie^-^< -^-L- exp{-<} . (48) 

For p2 approaching p^ one sees from Eq. (^ that C2 and C3 diverge, but cancellations occur. 
To take these into account we use Eq. (Bsl) to write 



C2e^-p + C3e^3.p = _ \ jA^' {P3(P3-R)e^--''-P2(P2-pi)e^^-n .^g. 

P2P3 iP2 - Pi ) iP3 -Pl) P2-P3 

and by recasting the expression in curly bracket as 

{iP3 - P2){P3 +P2- Pi)(f'^' + (P2 - P2Pi)(eP'"^ - e^'^'p)} (50) 
we obtain, with Eq. (^2|), 

Iri' ^40^ r / Q \ pP2Tp _ pPsTp ^ 

P2P3{P2 - Pl){P3 - Pl) l\2 / P2-P3 ) 

(51) 

which remains manifestly finite for p2 — P3. 

If all Pi are real, then it is evident that the expression in the curly brackets as well as the 
fraction in front are positive so that the overall expression is negative. Hence its omission leads 
in this case to an upper bound for Pq. 

If p2,3 are complex, ^2,3 = —a T ^b, then the last equation can be written as 



C2e^'2^p + cae^'^^p 



{a? + 62) (a2 + 62 + 2api + p\) 



I + 2p^ cos6Tp - 6sin6Tp + (a^ + api) g-'^^p (52) 
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From Eqs. (45) and ( [l6| ) for a and one obtains 

e-- < exp |-iVe -{^A+Ip,- r^} . (53) 

In Eq. (^) we bound the absolute value of the expression in the curly brackets by using 
|sin6rp/6| < Tp and | sin 6rp | , | cos 6Tp | < 1. With Eq. (|^) this leads to an upper bound of 
the form 

e-^*]e-^<= .r/' const [{a + /?rp)e~T"p] e'^'^ . 

The maximum of the expression in the square brackets is taken for Tp = Tm = - — ^ if Tm ^ 0, 
otherwise for 0. One can check graphically that Tm < for r]' < 1/3 and hence the maximum 
is taken at Tp = 0. An upper bound, valid for all A and 17, is then found graphically, with 
Ne = N^/i and for ry' < 1/3, as 

|c2eP2^p + C3eP3^p| < O.SlV e-^^/V _ (54) 
Together with the bound for cie^^^p in Eq. (^) this gives, for rj' < 1/3, 

^'o^(-p;|l)) < {rJl^ +0.81,'e-(^ -O^^je-^^ . (55) 

The expression in the curly brackets can be evaluated graphically, and this leads, for > 1 
and 77' < 1/3, to 

PoVp;I1))< (1+ e-^'^^^-l^^) . (56) 

In a similar way one can also obtain a lower bound for Pq, but it is quicker to use Jensen's 
inequality p^] which gives 



N=0 



With N = Ne Eq. (|2^) then yields 

(1 - = ei'^^i-"')^^ < . 

Since |ln(l - r/')| < ??' + i one finally obtains, with = N^iTp, |1)) and P^ = Po{Tp; |1)), 

g-(i+V/2)7Vd < (1 _ ^^./v' <p^<(i+ I r/')e-^<^ (57) 

3 

where the l.h.s. always holds and the r.h.s. for > 1 and < 77' < 1/3. 

Eq. ( |57| ) not only shows that ~ e~^'* is indeed an excellent approximation for small 
detector efficiencies, as already suggested in Section 3, but also gives strict error bounds. 



5 Conclusions 

A short laser pulse driving the 1-3 transition brings about a reduction of a state of the 1-2 
system in two ways as follows. If one or more fluorescence photons are emitted then the system 
ends up in |1), after a transient time for the decay of level 3. If no photons are emitted then the 
reduced (conditional) time development pushes the state towards |2). Its tiny |1) component 
can be estimated from above by 

|(ijV^)|V^^ . 
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Thus the reduction of \tp) to |2) is excellent if the parameters of the pulse are such that A^e > 5, 
i.e. if the average number of photons emitted by an atom initially in the ground state is 5 or 
more. If one measures with a detector which does not detect all photons (ry' < 1, e.g. solid angle 
less than in), then one has a problem for single systems. 

Detecting no fluorescence photons during a pulse may be due either to no emission at all or 
to missing the emitted photons. As a consequence a system with initial state for which no 
fluorescence photons are detected during a pulse is in a mixture of [1) and |2), namely 

K2|^)P \2){2\ + \m\'P,^{T,;\l)) . (58) 

Pq(tp; [1)) is the probability of detecting no fluorescence photons for a pulse, with the system 
initially in the ground state. We have expressed Pq in terms of the easily measurable quantity 
Nci{Tp; |1)), the average number of detected photons per system, with initial state is 
simply determined from the fluorescence rate of an ensemble of systems initially in the ground 
state. For small efficiency r]' we have shown that 

Jq — e 

and, more precisely, for r/' < 1/3 and A^d ^ 1; 

3 

Thus for increasing A'd the "non-reduced part" becomes very small, e.g. for = 3 it is already 
less than 5 %. 

Depending on the desired quality of the reduction one can adjust the parameters of the laser 
pulse in a simple way to give the required value of A'^d. 

If one considers not a single system but an ensemble of atoms the detector plays no role. After 
the pulse the state is reduced to a mixture of |1) and |2) regardless of whether all fluorescence 
photons are actually observed or not. This can also be seen by the use of Bloch equations 1^, p^ . 

If the 1-2 transition is driven by an additional interaction during the short laser pulse then 
complications arise and modifications to the projection-postulate result as obtained above have 
to be incorporated |^, |^. 
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Figure 1: V system with level 2 and auxiliary level 3 with Einstein coefficient A. $7 is the Rabi 
frequencies of the probe pulse of length Xp. 
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Figure 2: Typical form of fn'{p)- With increasing rj' the abcissa shifts downwards and passes 
for r/' = 1 through the point of inflection at pinfi- 
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